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ABSTRACT 



66 IRAS 23385+6053 is a Young Stellar Object with luminosity - 1.6 x 10 4 L 

at a kinematic distance of 4.9 kpc. This candidate precursor of an ultracompact 
Hn region is associated with a millimeter source detected at the JCMT but is 
undetected at centimeter wavelengths with the VLA. We observed this source 
with the OVRO millimeter array at 3.4 mm in the continuum, HCO + (l — > 0), 
H 13 CO + (l -> 0) and SiO(v=0, 2 1) line emission, and with CAM aboard ISO 
at 6.75/iin and 15/zm. The IRAS source is coincident with a 3.4 mm compact 
(r core ~ 0.048 pc) and massive (M ~ 370 M ) core, which is undetected at 15//m 
to a 3o" level of 6 mJy; this is compatible with the derived H2 column density 
of ~ 2 x 10 24 cm~ 2 and the estimated visual extinction Ay ~2000 mag. We find 
L S ubmm/L>boi ~ 3 x 10~ 3 and M eTO /M^. 3> 1, typical of Class objects. The source 
is also associated with a compact outflow characterized by a size r core , a 
dynamical timescale of iS 7 x 10 3 years, and a mass loss rate M £ 1O _3 M yr . 
The axis of the outflow is oriented nearly perpendicular to the plane of the sky, 
ruling out the possibility that the non-detection at 15 /im is the result of a 
geometric effect. All these properties suggest that IRAS 23385+6053 is the first 
& . example of a bona fide massive Class object. 
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1. Introduction 

The last few years have seen a rapid growth in observations aimed at identifying 
intermediate- and high-mass star forming sites in a wide range of evolutionary stages 
ranging from "Hot Cores" (Cesaroni et al. 119941 ), to ultracompact Hll (UCHll) regions 
(Wood & Churchwell |1989|) , to proto-Ae/Be stars (Hunter et al. |1998| ). The characterization 



of the earliest stages of high mass-star formation is more difficult than for low-mass objects, 
given their shorter evolutionary timescales. The likely candidates must be luminous, 
embedded in dense circumstellar environments, and non-associated with Hll regions (of. 
Habing & Israel pTO|) 



We have undertaken a systematic study to identify a sample of massive protostellar 
candidates, selecting sources from the IRAS-PSC2 catalogue and filtering them according to 
the above criteria. On the basis of H2O maser (Palla et al. |1991|) , ammonia lines (Molinari 
et al. |1996|) , centimeter and submillimeter continuum observations (Molinari et al. |1998a| ; 



1998b|) , we were finally left with about a dozen candidates which may be precursors of 



UCHll regions. Several of these sources have been recently observed with the IRAM-30m 
telescope (Brand et al. 1998, in preparation) and some with the Infrared Space Observatory 
(ISO), to characterize them over a wide range of wavelengths. We are now beginning a 
program of high spatial resolution millimeter wave observations of these objects. This 
Letter presents the results for the first source imaged, IRAS 23385+6053, at a kinematic 
distance [] of 4.9 kpc (using the Brand & Blitz (|1993|) galactic rotation curve). According to 



the definition (e.g. Andre |1996|) , this is an excellent candidate of a high-mass counterpart 



of Class objects. It is associated with an H 2 maser and with NH 3 line and submillimeter 
continuum emission, while no centimeter continuum emission was detected in our VLA 
observations down to a level of ~ 0.5 mJy/beam. 



2. Observations 

The field centered on IRAS 23385+6053 was observed at 88 GHz (continuum and 
lines) with the Owens Valley Radio Observatory (OVRO) millimeter wave array during 
January/February 1998. Two configurations of the six 10.4-m telescope antennas provided 
baselines in the range 15-240 m. Cryogenically cooled SIS receivers have typical average 
system temperature of ~ 350 K. The continuum observations used two 1 GHz wide 



1 In Molinari et al. (1996), we quoted a distance of 6.9 kpc derived from the observed galactic velocity 
field, which should however not be used for d>5 kpc in the 2 nd and 3 rd quadrants. 
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bands of an analog correlator. For the line observations, the digital correlator was 
split into several bands, enabling simultaneous observations of the HCO + (l — > 0) and 
H 13 CO + (l — > 0) rotational lines with ~ 20 kms -1 bandwidth and 0.84 kms" 1 resolution; 
the HCO + (l — > 0) line was also observed with a ~ 80 kms -1 bandwidth and 3.4 kms -1 
resolution. Additionally, data on the SiO(v=0, 2 — > 1) line were acquired in the image 
sideband of the HCO + (l — ► 0) line with the same bandwidths and resolutions. Gain 
and phase were calibrated by means of frequent observations of the quasar 0224+671. 
3C273 was used for passband calibration. The flux density scale was determined by 
observations of planets and the estimated uncertainty is less than 20%. The synthesized 
beamsize is 4 / .'l x 3'/5. Noise levels in the maps are ~ 0.7 mJy/beam, ~ 20 mJy/beam, 
and ~ 45 mJy/beam for the continuum and the low and high resolution spectral data, 
respectively. 

ISO (Kessler et al. |i"996l) observations with CAM (Cesarsky et al. |199^ ) were made 
in staring mode on June 18 1997. Diffraction-limited images (resolution 2". 2 and 5", 
respectively) were obtained in the two broad-band filters LW2 (A e //=6.75/im, AA=3.5/xm) 
and LW3 (A e / /= 15.0 jum, AA=6/im). The pixel-field-of-view is 3", and the absolute pointing 
accuracy of the satellite was presumably within ~ 3" at the time of the observations. 
More than 50 frames per filter were obtained with a 2s integration time per frame, and 
the expected 3a sensitivity limits were 0.33 and 0.45 mJy/arcsec 2 in LW2 and LW3, 
respectively. The data have been reduced using CAM Interactive Analysis (V2.0), and the 
absolute calibration is accurate within 20%. 



3. Results 

Figure |I|a shows the contour plot of the HCO + (l — > 0) integrated line intensity overlaid 
on the CAM-LW3 filter 15 /xm image. Figure |I]b shows an enlargement of the central 
region, with a contour plot of the 3.4 mm continuum. With the exception of HCO + (l — > 0) 
emission in one area to the NE, the HCO + line and the 3.4 mm continuum arise from 
the central region where the 15 fim continuum is faintest; HCO + also shows an extension 
towards S-SE. 

In the dark region at the center of the 15 /im image, flux densities of the order of 
~1 mJy/arcsec 2 are found. Similar values are also found in the southern portion of the 
CAM field, suggesting that such a constant level is due to faint foreground diffuse emission. 
The rms noise computed in the southern 20" of the image is ~0.15 mJy/arcsec 2 . Similar 
flux densities and rms levels are obtained for the 6.75 /zm image. We thus conclude that no 
15 fim source is detected at the position of the central 3.4 mm continuum peak down to a 
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3a level of ~ 0.45 mJy/arcsec 2 which, integrated over the deconvolved area of the 3.4 mm 
core, corresponds to ~ 6 mJy. 

The 3.4 mm continuum core is compact with slight elongations to the SE and N 
(Fig. [l]). The total integrated flux density is 19 mJy and the deconvolved source size 
is 4".5x3".6, corresponding to an average radius of 0.048 pc at a distance of 4.9 kpc. 
The continuum emission peaks at a(1950) = 23 h 38 m 31 s .4, 5(1950) = +60°53'50", and is 
probably due solely to dust thermal emission since the extrapolation at 3.4 mm of the 2 cm 
upper limit with an ionized wind model ~ z/ 6 (Panagia & Felli |1975|) yields less than 10% 
of the observed flux density. 

Fig. | shows the complete spectral energy distribution (SED) of IRAS 23385+6053. The 
OVRO flux density is consistent with those from JCMT (Molinari et al. 1998b ), indicating 



that no significant diffuse millimeter emission has been missed by the interferometer. On 
the other hand the IRAS flux densities include contributions of all the sources seen in the 
15 /iin image of Fig. |l|. Indeed, the IRAS 12 /mi and 25 /im points in Fig. [| (5.05 and 
17.8 Jy, respectively) are compatible with the CAM 6.75 and 15 /im flux integrated over 
all the bright emission visible on the CAM image (~4.8 and 9.3 Jy respectively, in an area 
of 6.8 xlO -8 sterad), while they lie more than 3 orders of magnitude above the CAM flux 
density limits for the core estimated above. Clearly, the 6.75 and 15 /im emission area, 
which is not directly related to the millimeter core, is responsible for the IRAS 12 and 
25/im flux densities; however, it is unlikely that it contributes for a significant fraction of 
the 60 /im and 100 /im IRAS flux densities because: 1) at millimeter wavelengths, the core 
is compact and no discrete or diffuse continuum emission is detected at either OVRO or the 
JCMT; 2) a power law extrapolation of the 6.75, 12, 15 and 25 /im flux densities to longer 
wavelengths, leads to flux densities at 60 and 100 /im < 10% of the observed values. We 
conclude that the millimeter core is the only source contributing significantly at 60 and 
100 /im. 

In Fig. ||the dotted line is a fit to the data, obtained adopting a spherical envelope model 
where density and temperature vary according to radial power laws with exponents -0.5 and 
-0.4 respectively. The density and temperature at the external radius, which is the observed 
i core ~ 0.048pc, are 7xl0 6 cm -3 and 40 K. The opacity is k u = Ki. 3mm (A(mm)/1.3)~ L9 , 
Ki. 3mm = 0.005 cm 2 g _1 (Preibisch et al. |1993| ). We compute the bolometric luminosity 
of the source by integrating the data with a power law interpolation scheme, excluding 
the IRAS 12 and 25 /im flux densities and using the CAM upper limit at 15 /im. At the 
adopted distance of 4.9 kpc, we find Lw ~ 1.6 x 10 4 L . Integrating the fitted density 
power law and adopting a gas to dust ratio of 100 by mass, we estimate a mass of ~ 370 M 
and a mean H2 column density of ~ 2 x 10 24 cm -2 for the core (corresponding to a visual 
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extinction Ay ~ 2000). We also computed the mass of the core from HCO + line emission. 
The optical depth in each channel within 3 km s -1 of Vlsr = ~ 51 km s" 1 was estimated 
from the HCO + /H 13 CO + ratio, assuming a 12 C/ 13 C= 87 abundance ratio using the galactic 
abundance gradient and the local-ISM value given in Wilson & Rood ( |1994| ) and a source 
galactocentric distance of ~ 11 kpc; we find that the core of the main isotope line is 
optically thick. We derive a lower limit (due to the optical depth) for the H 2 mass of 180 
M Q assuming T ex = 30 K, based on the peak HCO + (l — > 0) synthesized beam brightness 
temperature, and [HCO + ]/[H2]= 10 -9 . This value compares reasonably well with a virial 
mass of 160 M obtained from the line width (~ 3.5 kms -1 ) and HCO + (l — > 0) core radius 
(~ 0.06 pc). 

In Fig. [5] we present profiles of the SiO(v=0, 2 — > 1) and HCO + (l — > 0) lines at low 
spectral resolution and the HCO + (l — > 0) and H 13 CO + (l — > 0) lines at higher resolution, 
integrated over the ~ 20 arcsec 2 centered on the core peak. Both the SiO and the HCO + 
profiles show broad wings, while no wing emission is detected from H 13 CO + (l —> 0). Images 
of the HCO + and SiO emission integrated over the blue (from —63 to —53 km s _1 ) and the 
red (from —47 to —37 km s _1 ) wings are shown in Fig. ^. A compact outflow centered on 
the continuum source is detected; both the blue and the red lobes are visible in SiO(v=0, 
2 — >• 1), whereas only the blue lobe is clearly detected in HCO + (l — > 0), although a faint 
red lobe may be present. The lobes are barely resolved; their position relative to each other, 
i.e. their degree of overlap, suggests that the inclination of its axis can be no more than 
30° with respect to the line of sight (Cabrit Sz Bertout |1986j ; [19901) . This orientation of the 
outflow rules out any possibility that the large extinction derived toward IRAS 23385+6053 
and the non-detection at 15 /zm may result from geometrical effects (e.g. an edge-on 
disk) . Properties of the IRAS 23385+6053 outflow are listed in Table |l|. Column density, 
mass, momentum and energy are computed in each velocity channel and finally summed, 
while the dynamical timescale is averaged over the velocity channels. We estimated the 
column density of material assuming LTE and optically thin approximations in all the line 
wings channels within the velocity ranges listed previously, assuming the above determined 
excitation temperature of 30 K for both HCO + and SiO. For both molecules we adopt 
an abundance fraction of 10~ 9 relative to H 2 . The velocity measured in the lobes is a 
good estimate of the true gas velocity because the outflow axis is nearly perpendicular to 
the plane of the sky. This configuration prevents an accurate measure of lobe size, but a 
reasonable estimate can be made based on the non-detection at 15/zm. Since any outflow 
more extended than the core would have created a cavity through which mid-infrared 
radiation could have escaped, we infer that the extent of the flow lobes is less than the core 
size. The size of the lobes and the dynamical time scale of the outflow are therefore upper 
limits, while the mass loss rate and kinetic luminosity are lower limits. 
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The flow mechanical power and mass loss rate in Table [1] are high compared to values 
characterizing outflows powered by low-mass protostars (e.g. Fukui et al. |1993| ), but in good 
agreement with those for high-mass young stellar objects (Shepherd & Churchwell 1996| ). 
Nevertheless, the dynamical timescale is similar to values found for the youngest outflows 
around low-mass Class objects (Andre et al. |1993|). 



4. IRAS 23385+6053: A Massive Class Object 

IRAS 23385+6053 was not detected at 6 and 2 cm with the VLA B-array with 3cr 
sensitivity limits of ~0.5 and ~0.8 mJy/beam respectively (Molinari et al. |1998aj ), implying 



that any associated Hll region, if present, must be optically thick and very compact. In 
this latter case, from the VLA synthesized beamsize and assuming that the source has 
not been detected because of beam dilution, the 2 cm upper limit leads to an upper limit 
for the radius of an optically thick, spherical and homogeneous Hll region r# n < 100 AU 
(see also Molinari et al. |1998a|) . An Hll region originating from a BO ZAMS star (which 
has a luminosity comparable to the bolometric luminosity of IRAS 23385+6053) in a dense 
(hh 2 ~ 10 7 cm -3 ) environment takes only few years to expand to a radius of 100 AU (De 
Pree et al. |1995[ ). However, a modest accretion rate of 4 x 10~ 6 M /yr could squelch 



the expanding Hll region (see e.g. Walmsley 1995 ), thus explaining the properties of our 



source in terms of a heavily obscured B0-ZAMS with a residual accretion. Another possible 
explanation is that there is no Hll region because IRAS 23385+6053 has not yet reached 
the ZAMS. In this case we assume that the observed bolometric luminosity is due to 
accretion onto a protostellar core with accretion luminosity expressed as L acc = 3.14 x 10 4 
L Q (M*/M Q ) (R /R*) (M/10~ 3 M q yr _1 ). Since L acc is known, we can combine it with 
the mass-radius relation of Stahler, Palla & Salpeter (|1986l) , R* = 27.2 x R (M*/M ) - 27 
(M/10 _3 M Q yr -1 ) 41 , to derive the protostellar mass M* as a function of the accretion rate. 
For an accretion rate M ~ 10 -3 MQyr" 1 , which is expected if the outflow (see Table |Xp is 
driven by infall (e.g. Shu et al. 1988), we obtain M*=39 M , implying that the embedded 
IRAS source is a massive protostar. In either case, this source is evidently massive and 
extremely young. 

Based on the evidence presented here, we propose IRAS 23385+6053 as a prototype 
high-mass Class object it is undetected in the mid-IR (A < 15/im) and in the radio 
continuum; it has L su b mm /Lb i ^ 5 x 10~ 3 (where L su b mm is obtained by integrating 
longward of 350/im); it has a ratio of the envelope to stellar mass greater than 1 (see Andre 
et al. |1993|) , M env /M+ > 10. The very short dynamical timescale ( < 7000 yrs) of the 



associated compact molecular outflow is also compatible with values found for outflows 
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around low-mass Class objects. 
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CAM 15^m - OVRO HCO + (1-0) CAM 15/i.m - OVRO 3.4mm Cont. 
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Fig. 1.— a) CAM LW3 filter (15/an) image of IRAS 23385+6053, with superposed the 
contour plot of HCO + (l — > 0) line intensity integrated over a velocity interval of ~ 7 km s" 1 
centered at Vlsr — —51.0 km s _1 ; the CAM image has been resampled at the same 
pixel size used to reconstruct the OVRO maps. To avoid saturation at the peak emission, 
contours are spaced from 0.3 to 1.5 by 0.15, from 1.5 to 3.0 by 0.3, and from 3.5 to 5.6 
by 0.5 Jy km s _1 /beam. The dashed circle shows the OVRO primary beam HPBW. - b) 
Enlargement of the previous image around the central position with superposed the contour 
plot of the OVRO 3.4mm continuum emission from 2.0 to 9.0 by 0.7 mJy/beam. The circle 
shows the JCMT 19" aperture centered at the 1.1 mm peak position. The grey ellipse is the 
OVRO synthesized beam FWHP. 
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Fig. 2. — Spectral energy distribution of the IRAS 23385+6053 core. The dotted line shows 
the model fit described in the text. 



- 11 - 



i 1 r 



T 



i 1 r 



HCO + (1-0) 



1 - 



—3 



w 0.5 - 



i 1 r 



7xSiO (v=0,2-l) 



- 



J I L 



J I L 



J ,- - - - r 1 l 



-80 



-60 



-40 



-20 







1 1 1 1 1 1 1 1 1 
~_ HCO + (1-0) 

r-d" 


1 1 1 | 1 1 1 1 | 1 1 1 1 | 1 1 

i— ' .... 5xH 13 C0 + - 

| ^ 1 : 


" 1 1 1 1 1 1 1 1 1 


i i i 1 i i i i 1 i i i i 1 i i " 



-60 -55 -50 -45 

V lsr (km/s) 



Fig. 3. — OVRO spectra integrated over ~ 20 arcsec 2 area centered on the core. Top panel: 
low resolution, rms=0.02 Jy. Bottom panel: high resolution, rms=0.04 Jy. In all panels, the 
vertical scale is the flux density in Jansky. 
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Fig. 4. — Outflow maps of IRAS 23385+6053. Contour levels start at 0.3 and increment by 
0.12 Jy km s _1 for the two maps; the dotted contour corresponds to -0.3 Jy km s -1 . The 
cross represents the position of the peak in the 3.4 mm continuum map. The full ellipse is 
the OVRO beam 



Table 1. Outflow Parameters 



Parameter HCO+(l 0) SiO(v=0, 2 -> 1) 



N H2 (cm- 2 ) 


7.5xl0 22 


1.4xl0 23 


M H2 (M ) 


10.6 


19.6 


P = Y,iMiVi (M km s- 1 ) 


66.6 


159.6 


E=|EiM<«?(ergs) 


4.7xl0 45 


1.4xl0 46 


< T dyn > (yr) 

MCMoyr- 1 ) 


7300 


5600 


1.5xl0~ 3 


3.5xl0- 3 


P = P '/T dyn 


9xl0" 3 


2.8xl0~ 2 


L = E/Tdy n (L Q ) 


5.4 


21.6 



